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We give the details of the calculation of the spectral functions of the ID Hubbard model using 
the spin-charge factorized wave-function for several versions of the U — ► +oo limit. The spectral 
functions are expressed as a convolution of charge and spin dynamical correlation functions. A 
procedure to evaluate these correlation functions very accurately for large systems is developed, and 
analytical results are presented for the low energy region. These results are fully consistent with 
the conformal field theory. We also propose a direct method of extracting the exponents from the 
matrix elements in more general cases. 



I. INTRODUCTION 



After the recent photoemission experimentsba on quasi 
one-dimensional materials, the need of understanding the 
dynamical spectral functions of strongly correlated elec- 
tron systems has arised. While the low energy behav- 
ior is usually well described within the framework of the 
Luttinger liquid theory,0uu the experimentally relevant 
higher energies (« 100 meV) can be [Calculated for ex- 
ample by diagonalizing small clusterstj or by Quantum 
Monte-Carlo calculations. LI Unfortunately, both methods 
have limitations either given by the small size of the 
system or by statistical errors and use of analytic con- 
tinuation. Even for the Bethe-Ansatz solvable models, 
where the excitation spectra can be calculated, the prob- 
lematic part of calculating the matrix elements remains: 
The wave functions are required, and they are simply too 
complicated. There is, however, a special class of mod- 
els, where the evaluation of the matrix elements is made 
possible through a relatively simple factorized form of the 



wave- function, and sof. 
by Sorella and Pared! 
and by the authorsBE 



ie results were already published 
j for the insulating half-filled case 
away from half-filling. 
The dynamical, zero temperature one-particle spectral 
functions can be defined as the imaginary parts of the 
time ordered Green's function: 

A(k,u>) = — ImG(fc,cj), for id > fi; 



B(k,ui) = ImG(fc,w), for ui < fi. 

TT 

A{k, ui) is measured in angular resolved inverse pho- 
toemission experiments and can be calculated from the 
Lehmann representation: 

A(k, id) = | (/. N+l\alj0, N) j' 6(id-Ef +1 +E»), 



while B(k,id) is measured in the angular resolved pho- 
toemission experiments and is given by: 

B(k,id) = ^|{/,JV-l|a fc!CT |0,AT)| 2 5(a;-<+^f- 1 ). 

Here N is the number of electrons, / denotes the final 
states and a,k,a destroys an electron with momentum k 
and spin a. If the spectral functions are known, the time 
ordered Green's function can be obtained from 



G(k,id) 



, A(k,id') 

td — td' + i5 



" , B{k,d>) 

did 



ld — ld'—i5 



(1) 



The special models for which the matrix elements can 
be calculated are: 

i) The Hubbard model, defined as usual: 



H = —t 



U ^2 n ht n i,l> ( 2 ) 



in the limit U/t — > +oo; 

ii) The anisotropic t — J model 



H 



tJ 



-*£( 



liCT +h.c.) 



(3) 



i ct—x,y,z 



in the limit J a — > 0, where the usual projected 

operators. Actually, the Hubbard model in the large U 
limit can be mapped onto a strong coupling model usu- 
ally identified as the t — J modet-jpius three-site terms 
using a canonical transformation,EJ'll3 where J = 4t 2 /U 
is small; 

iii) An extension of the r£n— J model first proposed by 
Xiang and d'Ambrumenil,t2l defined by the Hamiltonian 



1 



i,cr 

+ E E ^(SrS^-^mni+^Vi,!, (4) 



where Vij = Y[j>=i(^ ~ n i+j') m the exchange part of 
the Hamiltonian ensures that two spins interact as long 
as there is no other spin between them. The motivation 
to study this model is that, unlike the infinite U Hub- 
bard model, there is a finite energy J associated with 
spin fluctuations, and this will give us useful indications 
about the finite U Hubbard model. 

From the models defined above, the Hubbard model 
is the most relevant one. It plays a central role as 
the generic model of strongly correlated electron sys- 
tems. Even though it is comparatively simple, it is 
very difficult to solve except for the one dipiensional 
case, where it is solvable by Bethe Ansatz.E-3 Unfor- 
tunately, the Bethe ansatz solution is not convenient 
for direct computation of spectral functions, therefore 
an alternative approach was needed. In the linait of 
small U one can use the renormalization grouptJ to 
show that the Hubbard model belongs to the. univer- 
sality class of the Tomonaga-Luitinger model,E3 usually 
referred to as Luttinger-liquid.Ej The Luttinger liquids 
are characterized by power-law decay of ccpselation func- 
tions, and nonexistence of quasiparticles.ta The under- 
lying conformal field theory can be used to relate the 
exponents to|.fmitc=sizc corrections of the energy and 
momentum.Mfi This gives consistent results not 
only with the renormalization group in the weak coupling 
regime,CJ but also with the special case of U/t — * +00, 
where the exponents of the static correlations, rpuld be 
obtained using a factorized wave function.ocjlHj 

Actually, the spin-charge factorized wave function also 
describes the excited states as welljHH and it can be used 
to calculate the dynamical spectral functions as well. The 
spectral functions obtained in this way are very educative 
and in some sense, unexpected. For example, it turns out 
that the spectrum contains remnants of bandstij cross- 
ing the Fermi energy at 3kp - the so called shadow 
bands. Also it gives information on the applicability of 
the power-law Luttinger liquid correlation functions The 
aim of this paper is not only to give the details of the 
calculation, that can be useful for other correlation func- 
tions, but also to present some new results on the low 
energy behavior of the charge and spin part (both for 
the isotropic Heisenberg and XY spin model). 

The paper is organized as follows: In Section [n] we re- 
view the factorized wave function and in Section |l| we 
show how the spectral functions can be given as a con- 
volution of spin and charge parts. Sections IV and |v| 
are devoted to the detailed analysis of the charge and 
spin parts. The relation to the results obtained from the 
finite-size corrections and conformal field theory is dis- 
cussed in Section VI. Finally, in Section VII we present 
our conclusions. 



II. THE FACTORIZED WAVE FUNCTION 

It has been shown,0-0 by using the Bethe ansatz so- 
lution, that the ground state wave function of the Hub- 
bard model in the U — > +00 limit can be constructed 
as a product of a spinless fermion wave function and 
a squeezed spin wave function \x}- This can. be alter- 
natively seen using perturbational argumentsEj and then 
extended to the t — J model in the J — > limit. More- 
over, the fv&ve function of the excited states are also 
factorizedl 



\N) 



K ({2})>®lx^(0,/fl)>. 



(5) 



The spinless fermion wave function \ip) describes the 
charges and is an eigen function of N noninteracting spin- 
less fermions on L sites with momenta 



kjL — 27t\Zj 



(6) 



where the Tj are integer quantum numbers and j 

1.2 Y. The charge part is not fully decoupled from 

the spin wave function \x), as the momentum Q = 
2ttJ/N {J = 0, 1, . . . , N - 1 ) of the spin wave func- 
tion imposes a twisted boundary condition on the spin- 
less fermion wave-function (each fermion hopping from 
site L — 1 to site will acquire a phase e % ® ) to ensure pe- 
riodic boundary conditions for the original problem. The 
energy of the charge part is 



E, 



N 



N 

-2t Y^cosfcj, 



(7) 



and the momentum reads Pj* = kj> or i using 

Eq. §: 



2tt 

T 



N 



N „ 2tt 



On the other hand, the spin wave functions |x) are 
characterized by the number of down spins Ni , the total 
momentum Q, and the quantum number /q within the 
subspace of momentum Q. They are eigenf unctions of 
the Heisenberg Hamiltonian 



N 



(9) 



i=l a—x,y,z 



with eigenenergies E s . J a depends on the actual charge 
wave function In the case of the U — * +00 Hubbard 
model, 



2f_ 1_ 
~U~N 



E {^\ n i n i+S- h \+5 n i h i-&Wl> ( 10 ) 



i,a=±i 
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where &t and b- are the operators of spinless fermions 

at site j. For the ground state \ip GS ) it reads J a = 
n(it 2 /U)[l - sin(27rn)/(27TO)], where n = N/L is the 
density. 

For the t — J model: 



(11) 



and for the ground state J a = J a n[l — sin 2 (irn)/ (7m) 2 ]. 
For the model of Xiang and d'Ambrumenil J a — nJ a 
and is independent of the charge part. The energy of 
the factorized wave function is then given as the sum of 
the charge and spin energies, with the assumption that 
the correct J is chosen. If £/ — > +00 or J — > 0, then 
the spectrum collapses and we can assume all the spin 
states degenerate, simplifying considerably some of the 
calculations to be presented later. 

Furthermore, we choose N to be of the form Al + 2 
(I integer), when the ground state is unique. Then 
in the ground-state the spinless fermion wave-function 
\4>l' ) 1S described by the quantum numbers Q = ir and 
{1} = {-N/2, N/2 - 2, N/2 - 1}, so that the distri- 
bution of the kj 's is symmetric around the origin and we 
choose the spin part as the ground-state of the Heisenbet-g 
model according to Ogata and Shiba's prescription.cil 
This choice of the spin wave function makes the differ- 
ence between the U — ► +00 and U = +00 (the so called 
t model) limits. 

The price we have to pay for such a simple wave func- 
tion is that the representation of real fermion operators 
a\ in the new basis becomes complicated. As a first 



r (l-n 3 -, s ) + ot 1 



step, we can write at a as at a = 

where at CT (1 — rij, s ) creates a fermion at an unoccupied 

site and the at ^rij^ adds a fermion at an already occu- 
pied site, thus creating a doubly occupied site, a means 
the spin state opposite to a. This latter process gives 
contributions to the spectral functions in the upper Hub- 
bard band, A UHB (/c, ui) which can be calculated in a sim- 
ilar way, but we will not address this issue in the present 
paper. 

Next, we define the operators Z ia and Z i acting on 
the spin part of the wave function: The Z\ a adds a spin 
a to the beginning of the spin wave function |xiv) if 
i = 0, or inserts a spin a after skipping the first i spins, 
and makes it N + 1 long, e.g.: z\ a \ ||) = \a Ti) and 
Z\<y\ TI) = I T u !)• The Z it(J is defined as the adjoint 
operator of Zj , i.e. it removes a spin from site i. 

Then, to create a fermion at the empty site j = 0, we 
need to create one spinless fermion with operator 6^ and 
to add a spin a to the spin wave function with operator 



aL(l-™o,ff) = Z\ b\. 



(12) 



creating a spinless fermions with b\ in the charge part, 
we have to consider the following two possibilities: cither 
the j = site is empty, and with a\ a we create a spin 

at the beginning of the spin wave function with Zq a ; or 
it is occupied, and we insert a spin between the first and 
second spin in \x) with Z\_ a . So we end up with 



(l-m, ff )= (1 - n )Zl a + n Z[ 



Obviously we choose the j = in further calculations for 
its simplicity. However, one can show that the final result 
does not depend on this special choice and the transla- 
tional invariance is preserved even for these complicated 
operators. 



III. SPECTRAL FUNCTIONS 

To use the factorized wave functions in the calculation 
of the spectral function it is more convenient to trans- 
fer the k dependence from the a k a operator to the final 
state: 



A(k,io) = J2 L \(f, N + 1 \4>JO,N) 



and 



y.5{u-E l ] +l + E£)5 kp N + i_ p r 



B{k,^) = Y J L\{f,N-l\a^%N) 



x 6 ( u - Eg + Ef - 1 ) S k p m _ P > 



where the momenta of the final states are Pf^ 1 ■ 

As we already pointed out, the addition of an electron 
to the ground state can result in a final state with or with- 
out a doubly occupied state. Correspondingly, the spec- 
tral function has contributions from the upper and lower 
Hubbard bands: A(k, oj) = A UHB (/c, w)-M LHB (/c, oA We 
will now consider A LliB (k,uj) only. From Eqs. (ra) and 
( |l2| ) we get the following convolution as a consequence of 
the wave function factorization: 



.4 



LHB 



(*,w)= C a {Q,<J)AQ{k,u>-<J), (13) 

Q,Uj' } tT 



The apparent simplicity is lost for a[ a . Then, apart from 



and similarly for _B(fc,w): 

B(k,w)= D a {Q,u')B Q {k,u-u'). (14) 

Q,uj' ,<j 

AQ{k,ui) and BQ(k,uj) depend on the spinless fermion 
wave function only: 
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A Q (k,u;)=Lj2\(^£({I})\bM£ S ) 
{i} 

x5(co - E^ 1 + E GS,cKp^-Pi B:C 
B Q (k,u;) = Lj2\(^. Q 1 ({I})\ b l<f 



x6(u-E» s>c + E»- l )6 kJ 



(15) 



and they are discussed in more detail in the next section 
(Sec. @). 

On the other hand, C a (Q,u}) and D <j (Q 1 oj) are deter- 
mined by the spin wave function only: 

2 



C„(Q,Lu)=J2\(XN + i(Q,f Q )\zl a 
fa 



x5(u> 



pN+l 
E f,s 



E GS,s)i 



D a (Q,u)=Y^\(xN-i(Q,fQ)\ZoMN S ) 

fa 

xSiu-E^ + E*- 1 ), (16) 

and are analyzed in Sec. [v| Although we do not present 
it here, a similar analysis can be made for A UHB (fc, cj). 

In Eqs. (13) and ( |l4|) the simple addition of the spin 
and charge energies is assumed. Strictly speaking, this 
is only valid for the U — > +00, J — > and the model of 
Xiang and d'Ambrumenil for any J. In the other cases 
the dependence of J on the charge wave function should 
be explicitly taken into account. Still, it is a reason- 
able approximation, as the important matrix elements 
will come from exciting a few particle-hole excitations 
only, which will give finite-size corrections to J in the 
thermodynamic limit. Furthermore, we are neglecting 
the t 2 /U corrections to the effective operators^ and to 
the wave functions. 



The momentum distribution function, = 



- I J. 



can be calculated from the spectral function as rik — 
J B{k,uj)doj, leading— to a similar expression as used by 
Pruschkc and Shiba:Ea 



n k 



= Y / B Q (k)D(Q), 



(17) 



(^^({^i^olV^Q' 3 ): where the two states have differ- 
ent boundary conditions. In the ground state Q' = 7r, 
but we will not specify Q' yet. To calculate these ma- 
trix elements, we need the following anti-commutation 
relation: 



1 j( k '-k)/2 r UQ'-Q)/2 sin ([Q' ~ Q}/ 2 ) 

L sin([fc' - k]/2) ' 



where k and k' are wave- vectors with phase shift Q/L 
and Q'/L, respectively, see Eq. (^|). For Q — ► Q' the anti- 
commutation relation is the usual one: {b\.,,b k } = 6k,k', 
while for Q =^ Q' the overall phase shift (Q—Q')/L due to 
momentum transfer Q — Q' to the spin degrees of freedopi 
gives rise to the Anderson's orthogonality catastrophe.Ea 
Then a typical overlap (0\b ktf . . . . . . b^JO), 

where |0) is the vacuum state, is given by the following 
determinant: 



{^fc^fci} {bk[,k 2 } 

{bk' 2 ,kx} { b k' 2 ,k 2 } 



{bk[,k N } 
{bk' 2 ,k N j 



{bk' kl } {bk'kz} ■■■ i b k',k N } 



Replacing the anticommutator, the determinant above 
becomes 



L -N e i(Q'-Q)N/2 -Q e ~i{k' ] -k ] )/2 gin 



. —l k', — fci • —l k\ — k 2 

sm — sm 2 — 



/V 



sm 2 n sm 2 n 



sm 



Q' -Q 
2 

— 1 k[-k N 
1 2 

— 1 k' 2 ~k N 



• ~1 kj.j—ki . — 1 k N — k 2 

sm N 2 sm w 2 



sm 



-1 k' N -k f 



where Bg(k) = j BQ(k,ui)du> and similarly D(Q) = 
fD(Q,uj)duj. 

The local spectral function A(u) = j^fc^f^) i s 
given by 



(18) 



where Aq{lo) = \ Sfc Aq(k, oj). Similar equation holds 
forSM. 



This determinant is very similar to the Cauchy deter- 
minant (there the elements are l/(fc — k') instead oL 
1/ sin(fc — k')) and it can be expressed as a product,Ej 
so for the overlap we get: 



±L -N e i(Q'-Q)N/2 sili 



N 



Q'-Q 



TT e -»(fc5-fcj)/2 



n 

j>i 



sm 



kj ki t r 

-2- It 



n 



V 



IV. ABOUT Aq(K,u) AND Bq(K,uj) 

To calculate Ag(fc, w) and BQ(k,u>) defined in 
Eq. (Ua), we need to evaluate matrix elements like 



where the sign + is for N = 1,4,5,8,9,.. and — for 
N = 2, 3, 6, 7,... 

Now we turn back to the AQ(k,u>). The matrix ele- 
ments in Eq. (|l5|) are 



4 



tL/,JV,GS\|2 



l' 



«^(m)i^i<o GS ; 



-22V sin 2JV 



Q'-Q 



j>l 



2 kj 



V — V 

2 «y 



it 



(19) 



where is a wave vector with phase shift Q' / L. Here we 
have used that 



JV , , , JV+1 



sin 



-1 h-q' 



T 2 ■ -2 Q' Q 

L sm — - — 



holds, independently of the actual quantum numbers {X} 
and {!'}. 

Similarly, for the matrix elements in E>Q(k,uj) we get: 



LM£({I})\b \t 



L,Q> 



= L- 2N+2 sm 2N - 2 



-2 Q'-Q 



n 



• 2 TT • 2 



j>i 



2 



We are now ready to calculate the spectral functions 
numerically. One has to generate the quantum numbers 
Xj, and evaluate the energy, momentum and the expres- 
sions above. 

From now on, we will consider Q' = it. 

First of all, it turns out that the following sum rules 
are satisfied for every Q: 

r * dk 



/ dk f 
— J dcoA Q (k,Lu) = 1- 

/ 2n J dujB Q{k^)=n. 



(20) 



In the absence of the Anderson orthogonality catas- 
trophe, when Q = Q' = n, the contribution to the spec- 
tral functions comes from one particle-hole excitations 
only, and the spectral functions are nothing but the fa- 
miliar 5{lo + 2icosfc). This is not true any more when 
we consider Q 7^ it. In that case we get contributions 
from many particle-hole excitations as well. The largest 
weight comes from the one particle-hole excitations, and 
increasing the number of excited holes, the additional 
weight decreases rapidly. Although from Eq. (|l^) we can 



calculate the matrix elements numerically for all the ex- 
citations of the final state, its application is limited to 
small system sizes (typically L < 30). It is due to the fact 
that the time required to generate all the possible states 
(quantum numbers X) is growing exponentially. There- 
fore, in some of the calculations we take into account up 
to three particle-hole excitations only. In Table. Q we give 
the total sum rule for small sizes in a calculation where 
we took into account up to one, two and three particle- 
hole excitations. We can see that the missing weight 
is really small in the approximation that includes up to 
three particle-hole excitations in the final state. So, if 
we restrict ourselves to a finite number of particle-hole 
excitations and introduce the function 



N/2 



7(1) = J] sin 2 ^(Z-Z') 

T' = -N/2 

N/2-1 . 

x J] sin- 2 i[l-l"]!L 

Z"=-JV/2 ^ 



2L 



(21) 



the calculation of the spectral weight becomes simple. 
The weight of the peak corresponding to a one particle- 
hole excitation can be given as: 



A Q (X p ,X h ) = 



1 



g(X h ) sm 2 ([X h -2?]£ ) 



-4, 



(0,0) 



(22) 



where we have removed the quantum number X h (hole) 
from and added X p (particle) to the set {X} of the ground- 
state of N ± 1 fermions, so that the momentum of the 



final state is P 
ergy is E? + 1 



N+l 
f 



= E, 



= kP - k h 
1 - 2t cos k p 



P, 



N+l 



and the en- 



2tcosk , where the 



GS 



JV+1 



GS — (N + 1)Q/L is the momentum of the ground 

state. Furthermore, Aq' ^ is the overlap between the N 
electron ground state with boundary condition n and the 
N + l electron ground state with boundary condition Q, 
and will be discussed later. 

Similarly, for the two particle- hole excitations we get: 



a t-rP tP Th Th\ 9(^1)9(^2) 
A Q{- L i,- L 2^ L i ^2) - 



sin 2 ([JC-^]f)sin 2 ([Xf-Xf]f; 



-A 



(0,0) 



5 (Jf) ff ^)sin 2 ([Xf-^]f)sin 2 ([Xf-^]f)sin 2 ([Jf-^]f)sin 2 ([X^^]f) « 
with energy and momentum P 



(23) 



JV+1 



h P 4- lc P h h h h -t- P N + 1 



GS 



E N+1 = E N +! 



2t (cos k p + cos &2 — cos k\ — cos k 



The corresponding equations for three or more particle- 
hole excitations are similar to those above, but since they 
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are long, we do not give them here. 

A typical plot of AQ(k,uj) is shown in Fig. |]. We 
choose Q = 7r/2, which is halfway between the symmet- 
ric Q — and the trivial Q = tt case. In the figure we 
can see the singularity near the Fermi energy, further- 
more the weights are distributed on a cosine- like band. 
To make it more clear, in Fig. || we show the support of 
Aq(H, u>) and the distribution of the weights. 



A. The weight of the lowest peak 

Now, what can we say about Aq'°\ the weight of the 
lowest peak? In the ground state the quantum numbers 
Xj and 2'j are densely packed, and from Eq. (|l9|) we get 



i (0,0) 



-,27V 



(0/2) 



L 2N 
N r 



N 

n 

3=1 



2 n J 



27V+l-2j 



><n 



. 2 (2j-l)n+Q . 2 {2j-l)n-( 
sm — sin 



2L 



2L 



j-N-l 



From this we can conclude that Aq' ^ is an even func- 
tion of Q and A^' Q ^ = 1. We are not able to give a closed 
formula for the sum. However, very useful information 
can be obtained by noticing that 



4(0,0) 



v sin 2 



.4 



(0,0) 



n 



2 2jn+Q ' 



"•Q-ir 3=1 

and in the thermodynamic limit 
(tt 2 - Q 2 ) 2 



4(0,0) 4 (o,o) 

^ L Q+27r^ L Q-27T 



(2Lsin7rn) 4 



1 



2L 



2tt 



cot7rn + 0(L 2 ) 



Here the Q is extended outside the Brillouin zone. Now it 
is straightforward to get the size- and filling-dependence 



of A 



(0,0). 
Q ■ 



(0,0) 



A 



where 



f(Q) 



(L sin7rn) Q Q 



a Q 777 cot 7771 
21j 



0(LT 



1 ( Q 



(24) 



(25) 



Eq. (f24j) is also valid for Bq°\ apart from the sign in 
the 1/L correction. 

The f(Q) is an even function of Q, /(tt) = 1, and it 
satisfies the second order recurrence equation 



f(Q + 2tt)/(Q - 2vr) (tt 2 - Q 2 ) 



2\2 



P(Q) 
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which can be reduced to 

/(Q + tt) r 2 (Q/2vr) 



/(Q-tt) r 2 (-Q/27r) 



and it follows that /(37r), /(57r) etc. are zero. In the 
interval from Q = to tt it can be approximated as 

O 2 O 4 
In /(Q) re -0.3047 + 0.3248-^ - 0.0201-^- 

with accuracy 0.0001. Furthermore m/(0) = -0.304637. 



B. Low energy behavior 

As we can see in Fig. ^ for low energies Ag(k,Lj) 
has so called towers of excitations centered at momenta 
k = (N + l)(Q + 2pn)/L, where p is an integer. The 
largest weights are for the peaks in the tower with p = 0, 
the next with p = — 1 (if Q > 0) or p = 1 (if Q < 0), and 
so on. The lowest excitation in tower p corresponds to a 
set of densely packed quantum numbers Xj shifted by p. 
From the definition of the momenta kj , this is equivalent 
to imposing a twist of wave- vector Q + 2pTt. Therefore we 
can introduce Q = Q + 2pn 1 where Q is not restricted to 
be in the Brillouin zone, but for p it has values out- 
side. We define Ag(k,uj) to describe the p-th tower, so 
that Aq(Jc, u>) has contributions from each of the towers: 
A Q (k,u)=J2 p AQ(k,u). 

Furthermore, we enumerate the peaks in a given tower 
with indices i and i' , so that the energy and momentum 
of the peaks are, from Eqs. @, (0) and (||): 

E i>v = E» s +e c + ^-u c (l + 9l)+^(i + i>), (26) 



2L 

Pi,i> = k Q + J + T^^' 



L 



(27) 



where we have neglected the 0(1/ L 2 ) finite-size correc- 
tions. Here e c = — 2icos7rn is the 'Fermi energy', u c = 
2isin7rn is the 'Fermi (charge) velocity' and kg = nQ 
is the 'Fermi momentum' of spinless fcrmions represent- 
ing the charges. By A^ 1 ■* we denote the weight of the 
peaks, and for convenience, we also introduce the rela- 
tive weights a { y' } = A i ^' ) /A^' 0) . The weight of the 
first few lowest-lying peaks can be calculated explicitly 
by Eqs. (^l|)-(|23|), as they are given by a finite number 
of particle-hole excitations. The degeneracy of each peak 
grows with i and i' . Here we assumed that the disper- 
sion relation is linear near the Fermi level with velocity 
u c . Clearly, this picture is valid for energies small com- 
pared to bandwidth. 

From Eq. ( p2[ ) we get the relative weights a^ 1 \ e.g. 

(1,0) • 
a ~ is given as: 
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1,0) 



sin 2 (^)sin 2 (^) 
S in 2 (f)sin 2 (^±2)- 



Introducing Wj — (Q/tt + j/) 2 /4, the relative weights in 
the thermodynamic limit simplify so that: 

a (0 ' 0) = 1 
a Q ■ 

(1.0) 



lowest lying peak in the tower determines the exponents 
of the correlation functions; ii) The weights in the towers 
are given by T-functionEJ 

The spectral function Aq(k,u>) in the thermodynamic 
limit is given by 



AQ{k,cj) = Y,A ( i' i ' ) 5(LO-E iti ,)6 k ,p it 



(32) 



(2,0) 

a~ 



2 2 



wi(w-i + w 3 ), 



and collecting everything together, Eqs. ( p4[ ) and (|3l]-[32|), 
for the low energy behavior of Ag(k, u>) we get 



(i,i) 

= W-xWx, 



A Q (k,uj)=Y, 



and also a 



Q 



Si' A 



f(Q)Q{Lo-u c \k-k Q \) 
47ru c sin(7m)r(/3Q + l)r(/3_Q + l) 



holds. Note that some peaks are 



degenerate and therefore they are a sum of more con- 
tributions. Now, it takes only one step to get the gen- 
eral formula which reads (including the finite-size correc- 
tions): 

( . _ (1 + /3q)(2 + /3q) ...(* + 
a Q ~ i\ 

(l + /3_q)(2 + /3_q)...(i' + / 3_Q) 
X V\ 



io — e, 



; + U c (k-k d )-\ f3 ' 



A-KUr sin nn 



lo-£ c 



-u c {k k{\) ^ 



4nu c sin(7m) 



(33) 



It is also worth mentioning the symmetry property 
oj) — A_q(— k,Lj). The whole calculation can be 
repeated for the spectral function Bq(k,ui): 



B Q (k,L0) 



f(Q)@(u c \k- 



^ 4TTu c sm(TTn)T(P Q + l)T(f3_Q + l) 



(i + i')n - (i-i')Q 



(28>< 



-UJ — u c (k+kp) 



A-KUr sin nn 



-Lo + u c (k+kq) 
AttUc sin irn 



(34) 



where 



2tt 



(29) 



It can also be expressed with the help of the T-function, 
since 

(1+/3q)(2+/3q) • • • (i+Pg) _ r(i+/?$+l) 



i\ r(z+i)r(/3 +i) 

The asymptotic expansion of the T-function gives 



(i + l/2 + /3 a /2)^, 



(30) 



r(* + i) 

which is a reasonable approximation apart from the i = 
peak. Then, it follows that aj~' 1 ' has a power law behav- 
ior: 

m + 1/2 + /3 Q /2)^ (i> + 1/2 + / 9_ /2)' 3 -« 



We should note, however, that these expressions are 
restricted for the weights far from the edges of the tow- 
ers, where the asymptotic expansion of the T-function, 
Eq. (j30|), is valid. This is especially true when Q — > ir, 
where the correct result is A 7r (k, u>) — S(uj — e c — u c [k — 
7m]). In other words, for the exponents close to — 1 
there can be a considerable deviation from the power 
law behavior, and the spectral weight accumulates along 
the edges of the towers. This behavior can be observed 
in Fig. [j], where the exponents are fi + q = —7/16 and 
/3_ Q = -15/16. 



1. Local spectral functions 



For the local (fc-averaged) spectral function Aq{u>) the 
weight of the j-th peak, denoted by A^ , is 



a~ — 



a(j) _ 1 A u 

A Q - T,2^ A Q 



T(P 6 + l)T(p_ 6 + 1) 



j'=0 



(31) 



The summation gives: 



Note that the exponent otq in Eq. (|2J) is also given by 

We can clearly see the manifestation of the underlying 
conformal field theory: i) The finite-size corrections to 
the energy and momentum [Eqs. (pq) and (p7|)] of the 



A 



U) 



1 r(l + qg+j) (0|Q) 
LT(l + aQ)T(l+j) Q 



tt n 2 — Q 2 0s 
1 + j 4- cot nn + 0(L ) 
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If we put it together with Eqs. ( p4| ) and fl2q), and ne- 
glect the 1/L corrections, the local spectral function in 
the L — > oo limit reads: 



1 



HQ) 



2nu c r(ag + 1) \27ru c sin7rn 



UJ — E<- 



(35) 



For Bq(u>) the ui — e c should be replaced by — ui + e c . We 
show Aq{lo) for some selected values of Q in Fig. ||. 



2. Momentum distribution function 

Here we try to make some statements about Bq(k) in 
Eq. (|l7|). A naive calculation in the low energy region is 
to sum up the weights near /ca 



E 

i=0 



5 



[»,-«+») 

Q 



if Z > ; 
if I < . 



Of course, one is aware that the summation includes high 

energies as well, where the equivalent for bf t , of Eq. ( ps| ) 
is not valid any more. However, the largest contribu- 
tions come from the low energy regions and the error is 
not very large. We do not want to get precise values, but 
rather some qualitative results. Neglecting the 0(1/ L) 
corrections, the sum gives for I > 0: 



D 



(0 



T(-a Q )T(l + l + f3_ Q ) 



For the J — > case the excitation spectrum of the 
spins collapse, and then we can use the local, ui inte- 
grated functions C a (Q) — '^Z u) C IJ (Q^uo) and D a (Q) = 
J2 W D a (Q, uj). They are related to the spirt-transfer func- 
tion ujji—yj^, defined by Ogata^and Shiba,EJ as was first 
noticed by Sorella and Parola.Q The spin transfer func- 
tion gives the amplitude of removing a spin a at site j' 
(here we choose j' = 0) and inserting it at site j , and 
can be given as 



wo-*a = (Xn \Pjd-i ■ ■■Pi,oS<,,s S \Xn ), 

where the operator = 2S.jSi+i + \ permutes the 

spins at sites i and i + Then C a (Q) and D a (Q) read 



C a (Q) 



N-l 



N + 1 



,i(Q+T)0'+l) 



LO . 



1=0 



N-2 



(37) 



3=0 



In particular, loq^o a = N a /N, and it follows that 
Eq C a (Q) = 1 and ^ Q D a (Q) = NJN. 

We are interested in these quantities for two particu- 
lar cases: the isotropic Heisenberg model because it is 
physically relevant, and the XF-model because it allows 
analytical calculations. We first consider the X F-model 
because the simplicity of that case makes it more conve- 
nient to introduce the basic ideas. 



and for I < the I and Q should be replaced by — I and 
—Q. Again, we can use the asymptotic expansion of the 
r-function to get 



. T(-a 



— sin(-7r/3 ± A) ' 



fcol 



^ V 2tt sir 



(36) 



where P_q for k > kg and /3q for k < kg should 
be taken in the argument of the sine. It is interesting 
that, although the exponent of the singularity oiq is the 
same for k > kg and k < kg, there is a strong asym- 
metry due to the prefactor (a, similar observation was 
made by Frahm and Korepino). In Fig. || this behav- 
ior is clearly observed. For Q — > ir the correct result of 
Bn(k) = 6(fc7r — k)Q(k 7T + k) is recovered. 



V. ABOUT THE SPIN PART 



To calculate C a (Q, uj) and D a (Q, to) given by Eqs. (|16[), 
we need to know the energies and wave-functions of the 
spin part. They can be calculated from the usual spin-i 
Heisenberg Hamiltonian, see Eq. (||), taking N and N±l 
sites (spins). 



A. XY model 

In this special case the spin problem can be mapped 
to noninteracting spinless fermions using the Wigner- 
Jordan transformation. It means that the eigenener- 
gies and wave functions are known, and we can calcu- 
late D a (Q,uj) and C CT (Q,w) analytically. We are fac- 
ing a similar problem - the orthogonality catastrophe - 
as when we calculated the Ag(w,fc), but now it comes 
from the overlaps between states with different number 
of sites. For convenience, we choose the spinless fermions 
to represent the a spins, so that the operator Z\ a (Z a ) 
only adds (removes) a site and does not change the num- 
ber of fermions, which we fix to be N s . Then we have 
to evaluate matrix elements like (xjv+i|-^o,o-IXjv S ) an d 
(XN-i\Z 0ta \x% s ), where in the |x° s ) the site is un- 
occupied and the fermions are on sites I — 1...N and 
from site I — 1 they hop to I — N skipping the / = site. 
For simplicity, we consider cases when the number of spin 
up and down fermions is odd (TV is even), so that we do 
not have to worry about extra phases arising from the 
Jordan- Wigner transformation. Then the momentum of 
the ground state |x^ S ) is -Pes = n - Let us denote by k' 
the momenta of fermions on a N ± 1 site lattice, quan- 
tized as k'j — 2irjJj/(N ± 1) and by k the momenta of 
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fermions on a TV site lattice, quantized as k = 2wJ r j/N, 
where Jj and J- are integers (j = 1, . . . ,N S ), and by / 
and p the operators of the spinless fermions. The energy 
and momentum of the state are: 



E = Jxy ^ cos kj , 

3=1 

3=1 



(38) 



(39) 



To calculate the matrix element in C (t (Q,uj), see 
Eq. (|l6|), we need the following anti-commutation rela- 
tion: 



{/fc"/fc} 



N N 



'N + i)Ti J-i 



ik I —ikl J 



VN(N , l=ll , =Q 

1 p j fc /2 sin(fcV2) 



V^V(JV + 1) sin([fc-fc']/2)' 



and the matrix element | (xn+i({J^'})\Zq Jx^) | 2 = 
\iO\fks, ■ ■ ■ fkjkjljk ■ ■ ■ fl Ng I0)l 2 ^ again given by a 
Cauchy determinant, which can be expressed as a prod- 
uct: 



2 n-j 



[n(n+ ir^n^f ny* 2 , 

3=1 3>i 



k' — k' 

| | sm 2 3 1 JJ i 



-2 &i % 



(40) 



3>* »>3 

Similarly, in the case of D a (Q 1 to), the anticommutator 



is 



{/£./*} = 



i 



y/N(N-l) 



,-ik'/2_ 



sin(fc/2) 



and the matrix element | {xn-i({J'})\Z 0jj \Xn S ) | i s 
equal to 



N„ 



3=1 

y — v 

i,3 



j>i 



sm 



-2 ^ ~ % 



(41) 



As soon as we have the product representation, it is 
straightforward to analyze the low energy behavior and 
also to obtain numerically D(Q, ui) and C(Q, ui) for larger 
system sizes. 



1. The low energy behavior 

The low energy spectra of D a (Q, u>) and C a (Q, u>) con- 
sist of towers centered at momenta Q r , a = 2m fi a , where 



r = 1/2,3/2, To analyze the low energy behavior in 

the tower labeled by r, we can proceed analogously to 
the charge part: the weights in the tower of excitations, 

Cr (t — — Cr,t7 Cr ,er and -Pr.cr = dr a Dr,cr , can be cal- 



culated from Eqs. fl40|) and (|4l|). The energy and momen- 
tum of the state can be calculated from Eqs. ( |3li| ) 
and ( j39|) and neglecting the 0(1/ N 2 ) corrections they 
read: 



E§f ] = E GS ±e a + ^u a ( 7 + a + 7 - CT + 2) 



(42) 



/':;■■"= Or,, ± ^ (7^ + (43) 



,(AT±1 

where 



A-(f^)'- 1 - 



(44) 



and the "Fermi energy" and the velocity of the spins are: 



£a = Jxy (^s cos it [i s - i sin n u, T ) 
u a = Jxy sin 717^, 



(45) 



and \x 9 =N S /N. 

The relative weights ; can be calculated from 
Eq. @, e.g.: 



7(0,1) 

1/2, o- 



,2 7r(l+Af+Af g ) -2 it(l+N+N s ) 

' bm 2(JV 2 -iV) 



JV-1 



jV-1 



2 7r(l+JV-JV,?+2JVJV g ) 



2(/V 2 -JV) 



7i 2 ((iV + iV ff + l)/2) 

i? 2 ((iV + JV ff -l)/2)' 



where 



Ar ff -i 



i?(Z) = JJ sin 



3=0 



7rZ 



7TJ 



N(N-1) NJ 1 



and the other cZr.o- are similar. In the thermodynamic 
limit, TV — > oo, the weight djyj CT simplifies to 

<il = (5 + f) 2 [i + o(inV£)]- 



(46) 



Neglecting the finite-size corrections, for general (i, i') 
and r we get: 

r(i + 7 - g + i) r(i' + 7 + g + i) 
r( 7 - CT + i)r(i + i)r( 7 + CT + i)r(i' + i)' 

where the exponents 7^. are defined in Eq. J44| ) and 
the weights again follows the prescription of the confor- 
mal theory, with strong logarithmic finite-size corrections 
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however. A similar analysis can be done for C <t (Q 1 uj). 
From the above and Eq. (jig) we obtain 



D a (Q,u>) « y^g(r, /y) [e CT - w + u a {Q - Q r>0 .) 

r 

x [e a - w - u a (Q - Q r ,a)] lr 
x6(e CT - a; + u fT \Q - Q r , a \) 



7 r 



(48) 



and 



C a (Q,uj) w [u3 — e a + u a (Q 



X [w - E CT - u CT (<9 - <5 r!CT )] 7r 

X9(w - £ CT - U CT |(3 - <3r,<r|) 



(49) 



where g(r, are numbers which can be determined nu- 
merically. 

We immediately see that the C a {Q) and D a (Q) are 
singular at Q — Q r ,a'- 

C a {Q),D a {Q) <x\Q - Q r ^° 
with exponent 

Vr,a = 7rtx + lr,a + 1 

and they are strongly asymmetric around Q r ,cn as we can 
conclude from the analog of Eq. ([36]) . 

For the non-magnetic case (/i CT = ji s = 1/2), the singu- 
larity is at Q r = ir/2 for all the towers, and the exponents 
of the main singularity (r = 1/2) are 7-^ 2 = —15/16 and 

7^2 = —7/16, furthermore 7^/2 = —3/8. 

B. Heisenberg model 

Although the Heisenberg model is solvable by Bethe- 
ansatz and in principle the wave functions are known, it 
is too involved to give the matrix elements of C a (Q,uj) 
and D„{Q,uj). The simplest alternative way-is exact di- 
agonalization of small clusters and DMRGlJ extended 
to dynamical properties. □ We have used both methods 
to calculate the weights for system sizes up to N = 24 
and N = 42, respectively. A typical distribution of the 
weights for C a {Q,uj) for zero magnetization is given in 
Fig. ||. There are several features to be observed: i) Due 
to selection rules, the nonzero matrix elements arc with 
the 5=1/2 final states only; ii) The weight is concen- 
trated along the lower edge of the excitation spectra in 
the interval 7r/2 > Q > 7r; iii) There are two, almost 
overlapping towers visible corresponding to r = 1/2 and 
r = —3/2. Our interpretation of the spectrum is that 
the weight mostly follows-jthe dispersion of the spinon of 
Faddeev and Takhtajan,E3 since the final states have an 
odd number of spins, thus there can be a single spinon 
in the spectrum and it has a cosine-like dispersion. It 
is also surprising that for C a (Q,uj) more than 97% and 



for D a (Q,uj) more than 99% of the total weight is found 
in this spinon branch. This behavior is suaailar to that 
discussed by Talstra, Strong and Anderson,t3 where they 
added two spins to the spin wave function. 

We can also try to analyze the low energy behavior 
from the conformal field theory point of view. Namely, 
from the Bethe-ansatz soU^kms.the finite-size corrections 
to the energy are knownE3 : L3o and they are also given 
by Eqs. ( j42] ) and ( |43| ) apart from ln(N)/N corrections, 
with 



7r,cr 



(50) 



For zero magnetization the velocity u a reads ttJ/2, the 
energy is e a = — J In 2 and £ = 1/V2, and the exponents 
are 7^, = — 1 and 7^, = —1/2, very close to the XY 
exponents (—15/16 and —7/16, respectively). For arbi- 
trary magnetization u ar —£ a and £ are to be calculated 
from integral equations.E3 

Also, we check if Eq. d47j ) is satisfied for the r = 1/2 
tower in Fig. [| Namely it tells us that c'- 1,0 -' = g?( 0,1 ) = 
1/2 and c^ 2,0 -* = d^ ' 2 ^ = 3/8, apart from finite-size cor- 
rections which we assumed to be of the same form as in 
the case of the XY model in Eq. (p6|). We believe that 
this method can also be used to determine exponents in 
a more general cases as well. 

Another interesting point is that the exponent J^m = 

— 1 already indicates that c' ' 1 ^ vanishes, in agreement 
with the selection rules. However, there is still some 
weight for c^ 0,2 -*, which comes from 5 = 1/2 bound states 
of spinons. We do not know the finite-size scaling of that 
weight, i.e. if it disappears in the thermodynamic limit 
or not. 

Now, if we recall that D a (Q) = J^u Ar(<2, ^), then 
it follows (see Eq. (|36|)) that the contribution to Dq for 
Q > 7r/2 is strongly suppressed, and we see essentially 
the contributions from the r = 3/2 tower. Since the 
contribution to C{Q 1 u>) and D(Q,u) come mostly from 
the lower edge of excitation spectrum, we can use the 
approximations 

C„{Q,u) = C ct {Q)5{w-s s -e q I 
D a (Q, u) = D a {Q)5{uj - e s + eq), 

where eq is the des Cloizeaux-Pearson dispersionB 



£ Q 



-J|sin(Q-7r/2) 



The C a (Q) and D a (Q) can be calculated numerically 
for small clusters (typically up to N — 26 with exact 
diagonalization and N = 70 with DMRG) for the non- 
magnetic case (see Refs. ||, |jj|). The (N + l)C a {Q) and 
(N — l)D a (Q) seems to have small finite-size effect, as 
follows from Eq. (|37j), and the singularity in the non- 
magnetic case is given hy rji/ 2 = —1/2, as already noticed 
by Sorella and Parolafl 
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We have also calculated C a (Q) and D a (Q) for the sys- 
tem with finite magnetization N±/N — 1/4 (see Fig. 0). 
There Qj = 37r/4, Qj = 7r/4 and the exponents are 
m/2A = -0.58 ± 0.03 and r] 1/2d = -0.25 ± 0.03. These 
exponents are consistent with £ = 0.87 ±0.02 and in sur- 
prisingly good agreement with the simple formula given 
by Frahm and Korepir£3 £ w 1 — /x;/2 valid in a large 
magnetic field. 



VI. THE GREEN'S FUNCTION AND THE 
COMPARISON WITH THE CONFORMAL FIELD 
THEORY 

The real space Green's function can be calculated from 
the spectral functions as 



G(x,t) = 



dk 



dtue 



iut — ikx 



A(k,u) 



and cd c ,d s are numbers. The quantum numbers D c , D s , 
AN C and AN S characterize the excitations and are re- 
lated to p and r as given in Table. III. The Z's are the 
elements of the so called dressed charge matrix. It can 
be calculated from Bethe Ansatz solution of the Hubbard 
model, and in the large U limit they read: 

Z cc = 1 Z cs = 
Z S c = Mi Z ss = £, 

where £ can be obtained solving an integral equation. For 
the non- magnetic case — 1/2 and £ = l/\/2. 

Then we are ready to identify the exponents: /3 ± q + 
1 = 2Af and 7^ + 1 = 2A^, and in this way we can 
directly see the validity of the CFT in the large- J7 limit. 

In case of the t — Jxy model no Bethe Ansatz result 
is known, but using the analogy with the isotropic case, 
the exponents are readily obtained using the substitution 
Z cc — > 1, Z cs — > 0, Z sc — > (i t and Z ss — > 1. 



for t > and A(k,u>) should be replaced by B(k, ui) for 
t < 0, as follows from Eq. (@). Then, from Eqs. @, @ 
and (|49|) it follows that: 



G{x,t> 0) 



E 



-iQ r xN/L 



1 

(a; - u s t)T^ +1 (a; + it c t)fr + 1 ' 



(51) 



where Q r was defined as Q r + 2irp, furthermore c p ^ r are 
numbers. The charge velocity u c is the same one as in 
Eq. (^), while the spin velocity is u s = u a Jn, where Mo- 
was defined in Eq. (|42|). The Green's function has sin- 
gularities at different momenta, depending on the actual 
quantum numbers p and r, see Tabic, [n] for details. 

On tho-other hand, according to the conformal field 
theory,cil'E3 a correlation functions {4>(x,t)(f)(0, 0)) reads: 



E 



c D c ,D a e 



-2i[D c fc T +(_D c +D s )fe i ]3; 



(x-u c t) 2A i {x + u c t) 2A < (x-u s t) 2A l (x+u s t) 2A ° ' 



D a ,D 

where the exponents 



VII. CONCLUSIONS 

To conclude, we have shown that for some special cases 
the spectral functions of the ID Hubbard can be cal- 
culated using the spin-charge factorized wave-function, 
which implies that the spectral functions are given as a 
convolution involving the charge and spin parts. Ana- 
lytical calculations are possible for the charge part and 
for the spin part in the case of the XY model. The low 
energy behavior turns out to be fully consistent with the 
predictions of the conformal field theory, i.e. the expo- 
nents are given by the finite-size corrections to the energy 
and momentum, and the weights are given by T-function. 
Based on this, we propose a new way to determine the 
exponents of the correlation functions. Furthermore, we 
argue that when the exponents of the correlation func- 
tions are close to integers, the Luttinger-liquid power-law 
behavior of the correlation functions should be taken with 
care, as it comes from the asymptotic expansion of the T 
function. 



, Z SS AN C - Z CS AN S 

2Af = Z CC D C + Z SC D S ± — 

1 2detZ 



2AJ = ( Z CS D C + Z SS D S ± ZccAN 2 S ~ et Z z ScANc ) <•*->> 
are related to the finite-size corrections: 
E E = ^u c (At + A") + ^u s (Af + A-) , (53) 



P - P Q = 2D c /c T + 2(D C + D s )k l 

+^(A+-A C 7 + A+-A S -). 



(54) 



* On leave from Research Institute for Solid State Physics, 
Budapest, Hungary. 
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TABLE I. Sum rule, Eq. (g0(), for Q = including one, two 
and three particle-hole excitations, N = L/2. 



L 


1 p-h 


1+2 p-h 




1+2+3 p-h 


4 


0.50000000 


0.50000000 




0.50000000 


12 


0.46477280 


0.49989083 




0.49999999 


20 


0.43436168 


0.49933463 




0.49999968 


28 


0.41165708 


0.49844924 




0.49999808 


36 


0.39388871 


0.49738700 




0.49999428 


44 


0.37941227 


0.49623473 




0.49998778 


52 


0.36725942 


0.49504054 




0.49997842 


60 


0.35682437 


0.49383182 




0.49996622 


TABLE II. The momenta for which the Greens function 


Gcr(x,t > 0) is singular. 
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TABLE HI. The correspondence between the Bethe 
Ansatz quantum numbers and p and r 



AN C 



p + r 
P 

—p — r 



—r 
r 
r 

—r 



12 



AQ(k,co) 



0.7 




k 



FIG. 1. A Q {k,uj) for Q = 48tt/97 (« 7r/2), N = 96 elec- 
trons on L = 192 sites. We can see the power-law singularity 
at k = 7r/4 and that the weight is accumulated along a co- 
sine-like band like structure. 
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~ n k Q -2jcn k Q 11 

k 

FIG. 2. Schematic plot of the support of AQ(k,ui) (above 
e c ) and Bq{— k,u>) (below e c ) for N/L = 1/3 and Q — tt/2. 
The dominant tower (p = 0) at k = fcQ and the sub-dominant 
tower (p = — 1) at k = kQ—2im are shown. The weight mostly 
follows the solid lines, and the shadowing represent the inten- 
sity. Although there are excitations above the dashed line for 
AQ(k,u>) as well, the weight associated with them is negligi- 
ble. The low energy part of AQ(k,ui) near k = kQ is enlarged 
on the insert, where the discrete states in the tower of exci- 
tations are shown. 



0.6 ■] Q=0 

Q=7i/2 --- 
0.5 Q=rc — 




0.5 1 t 1.5 1 2t 



CO 

FIG. 3. Aq(uj) for Q — 0, n/2 and n for quarter filling 
(L = 300, N = 150). For Q = the Van-Hove singularity is 
suppressed, and the weight is mainly near the Fermi energy. 
Q — 7r is equivalent to free-fermion case. The dotted line 
shows the low-energy approximation Eq. (|3^). 
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Q=0 

Q=447t/89 
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FIG. 5. The support and weights of C <y (Q,u>) for the 
N = 18 spin Heisenberg model. The symbols represent the 
excitations of the final states (19 spins), where the total spin 
is also indicated. The numbers near solid triangles give the 
weight of that particular state. Due to selection rule the ma- 
trix elements are zero with higher spin states denoted by open 
symbols. The dotted lines are a guide to the eyes and show 
the r = 1/2 and r — —3/2 towers. 



1 1 — — =i: . 1 1 1 1— -q^H, , 1 

-71 -27T/3 -71/3 tc/3 27T/3 71 

k 

FIG. 4. A Q (k) for Q = 0, 467r/91(fs tt/2) and 90tt/91(^ tt) 

(a) and B Q (-k) for Q = 0, 44tt/89(^ tt/2) and 88tt/89(^ tt) 

(b) for L = 270 and N = 90. The evolution of the weight 
and shape can be followed from the symmetric Q — case 
with the singularities at k — and k = ±27r/3 through the 
asymmetric Q — tv/2 case with singularities at k = -k/6 and 
— 7t/2 to the 'normal' distribution at Q — tt. 
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FIG. 6. The relative weights c (1 ' 0) , c (2 ' 0) , d ((U) and d (0 ' 2) 
as a function of the system size calculated by exact diago- 
nalization (squares and triangles) and by DMRG (crosses) 
for the r = 1/2 tower. The dashed line represents a fit to 
Oo+ai/iV+02 \og(N)/N form and it is reasonably close to the 
theoretical values 0.5 and 0.375 in the thermodynamic limit 
(a). The opposite sign of logarithmic corrections cancels if wc 
make the products [ c (i.o) d (o,i)]i/2 and [c (2,o) rf (o,2)]i/2 (b) 
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FIG. 7. C T (Q) and £> T (Q) (a), C|(Q) and D|(Q) (b) for fi- 
nite magnetization Nf/N = 3/4 with singularity at Q = 37r/4 
and Q = ?r/4, respectively. The solid symbols stands for 
D a (Q) and open for C CT (Q). 
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